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Abstract 

It is shown that step moving to meet solution flow can be unstable against lat- 
eral perturbations. The instability of long- wavelength perturbations occurs at 
values of the solution flow intensity less than some critical value depending on 
the step velocity. At given intensity of the solution flow, the instability comes 
at the step velocity exceeding a critical velocity. Decay of short-wavelength 
fluctuations is conditioned by the line tension of the step. The step moving 
along the solution flow is laterally stable at all values of the step velocity 
and the intensity of the solution flow. The overlapping diffusion field of the 
neighbour steps suppresses the lateral instability but it gives an instability 
of the step train against doubling of the period, i.e. neighbouring steps are 
attracted. The equidistant train moving to meet the solution flow is stable 
against the period variations. 
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I. INTRODUCTION 



The appearance of numerous investigations of step dynamics on the surface of growing 
crystals both experimental and theoretical is evidently stimulated by two factors. Firstly, 
development of technologies to obtain high quality crystal structures requires deep under- 
standing of elemental surface processes during the crystal growth. The control of the step 
dynamics plays a key role to provide the stable laminar growth, that gives high quality 
crystals. The experimental results had been mainly obtained by ex situ observations and by 
in situ very sensitive optical phase contrast or differential interference contrast microscopy. 
For the growth from solutions, observations of "live" mono-molecular steps on growing crys- 
tal were performed with the novel in situ optical technique Secondly, the advent of 
scanning tunneling and atomic force microscopy has open the door to investigate atomic 
scale structure of the crystal surfaces. In this measurements, the nanometer-scale precision 
can be achieved not only in the perpendicular to the surface direction as for the optical 
observations, but also along the surface. The direct observation during crystal growth from 
solution has been realized 0. 

Most studies have been devoted to the growth from vapour phases. Barton, Cabrera 
and Frank (BCF) in their seminal paper || who considered a step flow model, outlined 
this case. A vicinal surface consists of broad terraces separated by monoatomic steps. 
The transfer of the crystallizing matter includes adsorption upon the solid, diffusing across 
the terraces during some lifetime and further desorption back into the vapour phase or 
attachment to a step.Q BCF considered the steps as a perfect sink, when the attachment 
kinetics was neglected and the surface concentration of the substance had a local equilibrium 
value, taking into account the curvature correction via the Gibbs-Thomson effect. Taking 
into account of the kinetics of the adatoms attachment introduces, in general, two kinetic 
coefficients corresponding to adatom attachment from the lower and upper terraces [§]||. 



1 One dimensional diffusion along the step is also possible while the adatoms build into the solid 
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A linear stability analysis within quasi-static approximation |J predicted that the step can 
undergo a morphological instability at a step velocity of more than a critical value, only if 
the attachment rate from lower terrace exceeds the upper one. The instability mechanism 
is an analog of the Mullins-Sekerka instability [|7J. More general studies of the stability 
both the lateral and the longitudinal has been recently performed P-PI without using a 



quasi-static approximation. It was shown that the lateral instability can still take place for 
symmetric step kinetics because of an asymmetry in forward and backward direction of the 
step since it is in motion. 

The papers cited above are devoted to the linear stability analysis against small fluctu- 



ations. Scenario of further evolution of the step pattern was studied in Ref. JTT|]. This is 
evidence of presence of the spatiotemporal chaos in the step dynamics of the BCF model 
at the asymmetrical attachment kinetics. A competition between the thermal noise and the 
determinism was addressed in Ref. |T^| . 

Step bunching can be induced by an electromigration force during sublimation of vicinal 
surfaces, when the heating source is an electric current P~5|— |T5|| . Adatoms on the surface 



are assumed to carry electric charges and therefore their Brownian motion is modified by 
the electric field. The diffusion equation for adatom concentration c contains the term 
proportional FVc, where F is the electric force. The presence of the step bunching depends 
both on the temperature and on the direction of the heating current. In particular, reversing 
the current direction transforms a stable temperature region into an unstable one and vice 
versa. 

In contrast with the growth from gas phase, there are no convicting evidences on the 
mechanism of the substance transfer during growth from solutions. At the solution growth, 
the solute concentration is closer than the vapour density to the crystal density and the so- 
lution diffusivity is scarcely less than the surface diffusivity. Therefore, the volume diffusion 
in the solution with a direct incorporation into the step can be basic process transferring 



the crystallizing units flTB|. On the other hand, the direct incorporation can be suppressed 



by an interaction of the units with the solution and there is a possibility of the competition 



between the volume diffusion and the surface diffusion. In Ref. |17[], in situ observation of 
the dynamics of mono molecular growth step on the (1014) cleavage surface of calcite was 
presented and it was concluded that the surface diffusion does not control the calcite growth. 
However, the kinetic coefficient of these steps is very small ~ 2 x 10~ 7 cm/c, and the step 
motion is not controlled by any diffusion. The only attachment kinetics determines the step 
velocity. At the same time, there are reasons for the surface diffusion at growth of the KDP 



crystals []18|,|19|. We believe that investigations of step dynamics peculiarities at the volume 
diffusion can flow a light on this question. 

A formation of step bunches or a longitudinal instability forced by the influence of so- 
lution flow within the diffusion boundary layer was observed in Ref. [^] and has been 
considered theoretically in Ref. jH]] . The stepped vicinal was described as a continuous sur- 
face with kinematic waves of the step density. It was found that the step bunching depends 
on the angle of the riser and way of the heat transfer in the case of the melt growth. It can 



be developed only at the parallel motion of the straight steps and the solution. In Ref. [22 



the kinetics of the rectilinear macrosteps under diffusion or thermal interactions in stagnant 
media has been studied when the velocities of the steps are controlled by the supersaturation 
at the feet of the steps. The stability of the macrostep train depends on the angle of the 
step risers and way of the heat transfer in the case of the melt growth. 

The present study addresses the zig-zag instability of steps and the instability of the step 
train against doubling of the period during the crystal growth in the solution flow. 



II. GREEN'S FUNCTION FOR CONVECTIVE DIFFUSION 

Let us consider the diffusion nearby the vicinal crystal face being in a steady solution 
flow. We suppose that the height of the step shown in Fig. 1 is small in compare with the 
thickness of the diffusion layer and the characteristic curvature radius of the step, so the 
steps can be considered as a linear sinks for the diffusion field on the flat crystal surface. At 
large values of the Prandtl number Pr = vjD ^> 1, where v is the kinematic viscosity and 
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D stands for diffusivity, the thickness of the diffusion layer is much less than the thickness 
of the hydrodynamic boundary layer. The solution velocity can be developed in a series in 
perpendicular to the surface direction y. The solution velocity at the crystal surface vanishes 
and the velocity component parallel to the surface in the first order is 

v z = B{z)y. (1) 

The perpendicular component v y , being determined by equation of continuity 

dv y dv z 
dy dz 

has the form 

V 2 dB 

The component v y is in inverse proportion to a characteristic length of variation of the 
solution flow along the surface Lfi ow , i.e. in fact the size of the crystal face. Since the 
diffusion fields, varying in sufficient less scales, are of our interest, so the terms with v y can 
be omitted in diffusion equation and the parameter of the solution intensity B can be taken 
independent on z within the given part of the crystal surface. Then the diffusion equation 
for the solute concentration C takes the form 

where j(x, z) is the substance flux crystallized in this point of the surface. Let the solution 
velocity v z > and so B > 0, whereas the step velocity can be both positive and negative. 
After the Fourier transformation from t, x and z to u, p x and p z we obtain the equation and 
the boundary condition 

d 2 C _ _ dC 



(£o + £)C = 0; 



= (BD'ip^j, (5) 



where £ = (ip z A~ 2 ) 1/3 y and £ = (ip z A- 2 y 2/3 (icu / D + p\ + p 2 z ), A 2 = D/B. The solution of 
Eq. (||) is expressed through the Airy function vanished at y — > 00: 



C = 3 {BDH Pz )-^k\{i + eo)/Ai'(6). 



(6) 



Here at p z > 0, in order to satisfy the boundary condition at infinity we have to suppose 
Wz — Pz ex p{i' n '/2) and at p z < 0: ip z = \p z \ exp(— in/2). Taking these conditions into 
account, we obtain the Green's function of Eq. (H) in the form 



G(x,y,z,t) 



dt 



d,,, I dp z e iuJt+ipxX Re { e ipzZ ( — 

I \iPzj 



Ai 
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(7) 



Consider the diffusion field formed by the rectilinear steps moving as z = Vt at the 
constant solute flux j per unit of the step length. If Co is the solution concentration of the 
incoming flow, then 



- C(y, z-Vt) = -J- r d Pz Re e *"<- v *> ( — 
TlD Jo \ip z 



2\ V3 



X 



Ai 


{i Pz A- 2 


) 1/3 y + (zp z A- 2 )- 2 ^(pl-VD~Hp z ) 


Ai' 
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(8) 

Nearby the step, the concentration is logarithmical dependent on the distance to the step 
r = yfy 2 + C 2 ) C = z ~ Vt. The main contribution in the integral (§D goes from the large 
values of s = Ap z , when the ratio of the Airy function to its derivative tends as 



F(s) 



Ai 



{isf^y/A + (is)- 2 l*(s 2 - A 2 VD~Hs) 



_ s -2/3 e i7r/6-sy/A_ 



(9) 



Ai' [{ts)- 2 / 3 {s 2 - A 2 VD-Hs)] 

To use this asymptote as an approximation we have to correct the factor before exponent 
to provide the right value at s = 0: F = F(0) = -3 _1 / 3 r(l/3)/r(2/3), where T{z) is the 
gamma function. F(s) can be approximated by the function 



F{s) = -(s + s ) 

for which integrating in (H) gives 

J 



-2/3iir/6-sy/A 



So 



e W4| For 3/2 ; 



(10) 



Co - C = -^r(2/3)Re {* [2/3, 1; soA~\y - <)] } 



(11) 



Here ^(a, b; z) is the confluent hypergeometric function | 23| . In the limit r — > 0, we have 



C -C(r) = 4-ln^, a = 3- 1 / 2 
-kD r 



' (1/3) ,/2 e2^(i)^(2/3) ^ L90; (12) 



r(2/3)_ 

where ?/>(;?) = T'(z)/T(z) is the digamma function. The numerical integration gives cto ~ 
2.81. 

If the step height is small, li<C A, then the step velocity is controlled by the concentration 



(|i2|), where one has to put r « h. For the linear step kinetics, 

V = p[C e -C(h))/C e , (13) 

where f3 is the kinetic coefficient and C e is the concentration of the saturated solution, taking 
into account j = phV, we obtain 

V = ^ . (14) 

1 + pf3h{TiDC e )- 1 ln(a A//i) V ; 

Here a b = (C — C e )/C e is the bulk supersaturation and p is the crystal density. 

Let us investigate the asymptotes of the concentration field far from the step within 
quasi-static approximation. At y = and ( > in the limit ( — > 00, the expression 
(|TT|) is exact because, in this case, the main contribution of the integral (|Sp goes from the 
small values of p z , where the approximation (0) is exact. At ( ^> A, using asymptotic 



development of the confluent hypergeometric function [^] we arrive to the expression for 
the surface concentration 

To find the surface concentration at the side of flow incoming in the limit ( — > —00, we 
change the integration variable in (||) to r = ip z A: 

Ca . c(0 = ± lm f dTT -^M-^l (16) 

At ( < 0, the path of the integration in the r-plain shown in Fig. 2 can be reduced to the 
positive part of the real axis. The integrable function has poles only on this new path, where 



the function is real. So the imaginary part of the integral is the half-sum of residuals in 
the poles. If \£\ ^> A, the main contribution goes from the closest to zero pole r = (a^) 3 / 4 , 
where a[ ~ 1.01879 is the least root of the equation Ai'(-a') = 0. The result is 



Co - C(C) = 



1-3/2 



cxp 



C<o. 



(17) 



The concentration field relaxes slow along the solution flow and fast in the opposite direc- 
tion. Analogous calculation for the concentration far from the crystal surface leads to the 
expression 



C - C{y) 



3j 



(18) 



W 2 D (a' 1 ) 25 / 16 Ai(-c 

Note that the quantity A can be interpreted as an effective thickness of the diffusion layer 
in this case. 

At the large step velosities, V ~ .DA" 4 / 3 ^ 1 / 3 , the quasi-static approach is not correct. In 
the limits V ^> Z)A~ 4//3 £ 1//3 and |£| — > oo, the diffusion field reads as 

where w = VA/D characterizes the degree of the quasi-stationarity of the diffusion. Like 
the quasi-static case, the zone impoverished by the step sink is drifted by the solution flow. 
The results obtained in this Section will be used for analysis of the step train stability. 



C -C(()=jx 



(19) 



III. STABILITY OF ONE ISOLATED STEP 

A step can be treated as isolated if the diffusion length A is small in compare with the 
step spacing. Certainly, the isolated step can contain several elementary crystal layers. We 
shall assume that scale of internal structure, such as height and width of this composite step 
is the least of sufficient length parameters. For the isolated step, the scale plays a role of 
effective step height h which determines the concentration (|1^) controlled the step motion. 

Near by the crystal surface, the concentration up the solution stream is higher than the 
concentration down the stream. The isolines of the concentration field are schematically 
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shown in Fig. 3. It is clear that a negative local deviation of the step from a straight, 
Sz(x) < 0, being directed upstream, finds the more concentrated solute and the opposite 
deviation meets the impoverished solution. Further evolution of the perturbation depends 
on the relation between directions of the step motion and the solution flow. For the step 
moving along the flow, being at the less concentration the leading part of the step becomes 
slower. At that time, finding the larger concentration the dropping behind part overtakes 
the rest step. Therefore, the lateral perturbations of the step moving along the flow always 
decay. At the case of the opposite motion, the velocity increment of the perturbed part is 
having the same sign with the variation of the z-coordinate of the step. The presence of this 
positive back-coupling can produce a spontaneous growth of the perturbation, i.e. to the 
instability of the straight step shape. The line strain of the step and forcing the solution 
flow stabilize the straight shape. 

Thus, the flow of the solution gives rise an analogous effect to the Schwoebel's barrier 
HI at the surface diffusion transport. In the case of the volume diffusion, the substance 
flux is larger from the side of the solution incoming than from the opposite side. At the 
surface diffusion, the fluxes arriving from the opposite terraces are differ due to a difference 
of the energy barrier of the adatom attachment from the upper and the lower terraces. The 
feature of the volume diffusion is in a possibility to govern the quantity of the attachment 
asymmetry by the direction and the intensity of the solution flow. 

Consider the kinetics of small fluctuations, z = Vt + 5z(x, t), from the straight step. The 
unperturbed step velocity V is determined by (p~3|) . Taking into account the Gibbs-Thomson 
effect, the kinematics equation takes the form 

r i —3/2 

z = ±/3a(z,t) +[3\z" [1 + {z') 2 \ , (20) 
where p c (a) = A/cr is the critical nucleus radius on the crystal surface at the supersaturation 



a near by the given point of the surface and A is the capillary length. The upper sign ( |20| ) 
corresponds to the parallel motion of the step and the flow and the lower sign corresponds 
to the opposite motion. Assuming 5z(x, t) = e(t) cos(qx), let us calculate the concentration 



near by the step at eq <C 1. The substance flux per the unit length is j = ph [V + e cos(gx)]. 
Using the Green's function (|7|), we obtain 



a(x, t) = a b - (ttU)- l V \n(a /H) - eI (H, Q) + eVA~ l I{Q) U~ l cos(gx) 



(21) 



_1 p dr Ai[g(^)V 3 + (zr)- 2 / 3 (r 2 + g 2 
1 {H,Q) ^Ke^ ^ 1/3 Ai'[(ir)" 2 / 3 (r 2 + g 2 )] 



1 r 00 
/(g) = - Re / dr 

7T JO 
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Ai 
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Ai' 


[-(.r)V 3 


Ai' [(<r) 



( zr )-2/3 (r 2 + g 2 



(22) 



(23) 



where U = DC e /(ph), H = h/A and Q = qA. Taking into account (^T|) reduces the 
kinematic equation ( p0|) to 



/3 ^^f/- 1 / - AA- X g 2 

e = v e, v — — 
A 



(24) 



1 + PU-Hq 

where v are the amplification rates for the parallel and the opposite motion respectively. At 
positive v, the perturbation amplitude e is exponential growing in time and negative values 
imply decaying. The function Iq(H,Q) and I(Q) are always nonnegative. The dependences 
Io on H at several values Q are shown in Fig. 4. At H 1, they can be written as 



I (H,Q)=7f- 1 \n[a(Q)/H}. 



(25) 



The dependence a~ l {Q) is shown in Fig. 5, where a(0) = a , a(oo) = 2 exp(— 7) and 
7 ~ 0.5772 is the Euler number. The dependence I{Q) are shown in Fig. 6. At Q <C 1, 
we have I(Q) ~ 0.19g 2 and /(oo) w 0.128. With a relative error below 1% the functions Jo 
and / can be approximated by the expressions 



I Q (H, Q) = (27T)- 1 In [(0.77 + 0.12Q 2 )/H 
I(Q) = 0.19g 2 /(l + 1.58g 2 ). 



(26) 
(27) 



The expression (|24] ) shows that the step moving along the stream is stable for the lateral 
disturbances at all wavelengths. At the opposite motion, if I"(0) > UX/(VA), then the step 
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is unstable for the long wave disturbances with q < q cr , where q cr ps A 1 [0.12AV A /(A?7) — 
0.63] 1/2 . The instability condition is given by 

Vh < C e A£>V2- [ 8j 

Thus, the lateral instability occurs at low intensity of the solution flow or /and at sufficient 
high values of the step velocity and the step height. The instability zone, < q < q cr , of the 
disturbance wave vectors is shown by shading in Fig. 7. The thick line, q cr , starts at origin 
because the step is stable without flow in the frame of quasi-statics. The minimal wavelength 
of the unstable disturbances l mm = 82.7(U/V)X is determined by the capillary length times 
the ratio of the parameter of the diffusion velocity to the step velocity. Remember that the 
quasi-static approximation is correct at VA/D <C 1. 

At low intensities of the solution flow, B — > 0, the characterictic size of the impoverished 
zone A — > oo, therefore the instability vanishes in this limit. At high flow intensities, the 
concentration near by the step is close to the bulk concentration at all points of the step 
that explains for the stability at large B. Thus the instability vanishes at q = and q = q cr 
and the amplification rate of the instability has a maximum within the interval < go < Qcr- 
If at least one of the inequalities g A <C 1 and (3/(2irU) <C 1 is fulfilled then 

XU \ 1/2 



g A = 0.80 



'0.19AVV /2 J V2 0.12PV 



/0.19AKV 
\ XU ) 



- 1 



^max A.U 



0.19AV; 

The instability requires the step track length L = V/u max to manifest itself. 



- 1 



(29) 



The single step approximation addressed in this Section becomes inapplicable when the 
step spacing d is comparable to the diffusion length A. If a screw dislocation, having the net 
number m of unit steps in the normal component of the Burgers vector, generates a step 
train, each step consisting of n of unit steps, then 

d w 20np c /m. (30) 

At d < A, the single step instability can be realized by a sharp dropping of the solution 
temperature from a value at which d* ^> A in the steady state. The spacing d* will have 
been conserved while new generated steps have not arrived. 
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IV. LATERAL STABILITY OF VICINAL FACE 



To understand the influence of the diffusion interaction of neighbour steps it is appro- 
priate to consider the inverse limit to the single step case. At d <C A, the crystal relief 
is described by a height function h(x,z). Volume conservation implies that a continuity 
equation satisfies: 

§ + V(PV) = 0, P = vft , V = (31) 

where |P| is the surface slope to a singular face and the direction P is parallel to the step 
velocity V. The normal growth rate of the crystal face R = PV and the flux per unit 
area is J = pR. Using the Green's function (|7|), for the small fluctuations from the plain 
P = (0, P ) of the form 

P = P + P exp(iut + iqx + kz) (32) 

R = R + Rexp(iut + iqx + kz), R = PoV , (33) 

where Ro and V = (0, Vb) are unperturbed normal growth rate and step velocity, we obtain 
the surface concentration in the form 

C = Cb — pAD~ 1 T(u, q, K)Rexp(iut + iqx + kz), (34) 

(35) 



Ai 


/A- 2 *) 


-2/3 


[iujD 


l + q 2 - k 2 ) 


Ai' 


\A- 2 k) 


-2/3 


(iujD- 


- l + q 2 -K 2 ) 



These perturbations are the quite lateral only at k — 0. The stability against the quite 
longitudinal perturbations, q = 0, has been considered in Ref. ||21|| . Using the kinematic 
equation (|3l|) , we obtain the linear stability equation 

iuj + Vqk + Pc R q 2 - iT(u, q, K)ojpP Q (3A(DC e y l = 0. (36) 

In the limit k" k', where k' and n" are the real and imaginary parts of k = k' + in", the 
linear dispersion relations are given by 
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K " = -u/Vo = A _1 0, k' = PopC^pe Im(T) - P p c q 2 . (37) 

The amplification or decay lengths for the step motion up and down the solution flow (upper 
and lower sign respectively) are given by 

r L± r 1 -J p ° l «V/'iJrV3 a-[W- 2/3 W 2 + " 2 t^)] ) \P M 

y L ) ~ ± -c^d b lm \' s? ppTw + p i »«)] / " "~ A^— ' ( ) 

where w = |Vo|A/ZJ takes into account the influence of the step motion on the diffusion field. 
Positive and negative signs of L imply the instability and the stability respectively. For quite 



longitudinal displacements of steps from steady state positions, Q = 0, the expression (p8|) 
is a particular case of the results of Ref . |2l| . 

Quite lateral fluctuations, 9 = 0, are always decay, L 1 * 1 = — \P \p c q 2 < 0, the decay 
rate being determined only the step line tension. The explanation lies in the variation of the 
concentration along curved step being compensated with the diffusion influence of neighbour 
steps. 

It can be numerically shown that at the opposite motion of the steps and the solution 
flow, at all values of the parameters, the length L~ < 0, hence the strong diffusion interaction 
completely suppresses the lateral instability. In the case of parallel motion, there is bunching 
with simultaneous amplification of the perturbation along the steps for only sufficient long- 
wave perturbations: 6 < 9 cr (w,Q). The presence of the lateral perturbation, Q / 0, 
decreases the instability of the density waves of steps i.e. straight steps are more unstable. 
Note that the diffusion nonstationarity has a stabilizing character: the instability drops with 
increasing w. 

Thus, the diffusion interaction suppresses the lateral instability at A > Ad, where A ~ 1. 

For the dislocation steady growth (|30|) , it confines the unstable zone of supersaturation from 

above. As a whole, the unstable zone reads as 

5.3XU n fB\ 1/2 20A\n f B\ l/2 , . 

<°<— — hd , (39) 



fin \Dj m \D, 

where Uq = DC e / (ph ) and h is the height of unit step. This zone exists only at n 2 /m > 
0.26U /(A/3). 
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For an example, evaluate the instability of step on the (101) face of the KDP crystal. 
The value of the solution flow intensity is given by 

d 2 ip 



B 



(40) 

y=o 



dy 

where ip is the stream function. B can approximately be evaluated as B ~ U so i/S, where U so i 
characterizes the flow velocity far from the surface and S is the thickness of the hydrodynamic 
boundary layer at this point. For concreteness, let us take a tangential flow of solution. Other 
useful plane mass transfer models have been described in Ref. [24fl. For the flow around a 
semi-infinite plate z > 0, at an incoming velocity U so i, we have B = 0.33(Ug Ol / 'vz) 1 ! 2 . 
For U so i — 10 cm/c and z = 1 cm we obtain B = 100 c" 1 , and the instability zone is given by 
0.3/n < a < 0.003n/m, it existing if n 2 /m > 100. Taking into account an anisotropy of the 



step kinetic coefficient and the isotropy of the critical nucleus [IS] weakens this restriction, 



because the step spacing in the direction of the largest kinetic coefficient exceeds 20p c p6 



At the case of m — 1, the unit steps are stable against the lateral fluctuation but the steps 
containing ten or more units are unstable. At a = 0.05, for the step of 20 unit layers, the 
amplification length L~ pa 0.15 cm that can be observed. 



V. INSTABILITY OF STEP TRAIN AT WEAK INTERACTION 

The diffusion interaction between steps of a train is weak when the step spacing is 
considerably greater than the diffusion length: d ^> A. If the steps are moving to meet 
the solution flow then the condition of the lateral instability is given by ( p8|) with a weak 
stabilizing influence of the interaction. A sufficient effect of the weak coupling occurs at 
the parallel motion of the step and the solution. In this case, the lateral instability does 
not develop, so we will assume that the steps are straights with the coordinates z n . Since 
the diffusion field going from a step to meet the solution flow decays faster than along the 
flow, we will take into account only influence of the closest step from the side of solution 
incoming. The velocity of nth step is given by 
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2/3 



v, 



n 



V 



2nDC e \z n -z n -i) 



Jn—ly 



(41) 



where the unperturbed step velocity V is determined by (|Hj). The perturbation doubling the 



train period, z n = nd + (— l) n eexp(z/t), has a maximal amplification rate. The amplification 
length for parallel motion and the decay length for antiparallel are respectively given by 



The magnitude of x — {P/^U) hx(ao/H) is a criterion of the regime of the step motion. In 
the limit of the small values, the step velocity is controlled by the crystallization kinetics, 
i.e. the kinetic regime takes place, and at the limit x ^ 1 the diffusion regime occurs. Note 
that even near by the kinetic regime, when the step velocity practically does not depend on 
the flow velocity, the instability nevertheless exists though the amplification rate is small 
with respect to this parameter. 

The reason of the instability of the train period at parallel motion with respect to the flow 
is quite clear because if a step displaces forward then the step velocity takes an increment 
due to weakening the diffusion overlap of the previous step, and at the same time, the next 
step becomes slowing down, hence the perturbation will grow. This effect is an analogous 
to the current-induced step bunching during sublimation [p^ -[T5|] for the step up direction of 
the force if the surface diffusion length \2d is small compared to the step spacing d. In this 
case, the induced migration is directed along the step motion. However, the migration force 
can also induce the bunching for the step down direction if there is the opposite correlation 
between the surface diffusion length and the step spacing. This correlation depends on the 
temperature. Increasing the temperature, the surface diffusion length decreases exponen- 
tially and may become smaller then d. At the volume diffusion growth and dissolution, the 
bunching can takes place only for the parallel motion of the steps and the solution flow at 
any relation between the diffusion length A and d. At the surface diffusion, if the d is small 
compared to A 2 d, the most adatoms departed from a step reach the adjacent step to which 
the migration is directed. The adatom concentration grows along the force because they can 




(42) 
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only incorporate into the step or evaporate. It is precisely this fact that gives the instability 
for the opposite direction of the step motion and the adatom migration. On the contrary, at 
the volume diffusion, the solute reached the following step can move further with the flow 
and the opposite motion is always stable. 

At large step velocities laying out of the frame of the quasi-static approximation, the 
diffusion field of step has the power-type asymptotic along the solution flow and the expo- 
nential decay in the opposite direction fll9|) . Hence in this case, the step train is also stable 
against the doubling of the period during motion in the opposite direction with the flow and 
it is unstable during the parallel motion. 

VI. RESULTS AND CONCLUSIONS 

The linear stability analysis has shown that, at the opposite motion of step and solution 
flow, the lateral instability is possible, and, at the parallel motion, the instability against 
perturbation in the distances between steps can occur. The overlap of the step diffusion 
fields suppresses the lateral instability and gives rise the longitudinal instability of step train. 
Intensification of the solution flow and decrease of the step height stabilize the equidistant 
train of straight steps. Though the instability has a diffusion nature, it can be even noticeable 
near by the kinetic regime of the step motion. While the quasi-static approximation is 
correct, the lateral instability becomes stronger with an increase in the step velocity. Out 
of the quasi-static frame, the effect of step motion plays a stabilizing role. 

The evaluations have shown that the instability considered here can be observed under 
the real conditions of solution growth of crystals like KDP. As a rule at solution growth, 
meandering unit steps and presence of the straight macrosteps allow to grow high quality 
crystals. A danger of defect formation like solution inclusions does go from the presence of 
meandered macrosteps. The concave segments of such steps further a development of over- 
hanging upper step layers, and then the formed cavity can capture solution into the crystal. 
During the high-rate growth of large KDP crystals, the crystal surface with macrosteps is 
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commonplace, so providing with uniformity of the diffusion layer thickness p7| , |28[| does not 
guarantee a stable growth without the inclusion formation. Increasing the convection until 
the kinetic growth regime can be also insufficient condition to obtain high quality crystals. 
It is important to note that fast reversal of the flow direction can eliminate the instability 
both lateral and longitudinal. If the interruption time between equal intervals of the op- 
posite solution flow is small in compare with the time of the instability development 
then we find from (^4j) that the lateral shape perturbations decays with the average decay 
constant 

V=-(3Xq 2 {l + (3U- 1 h)- 1 . (43) 
A model of the diffusion layer relaxation when varying the flow velocity was considered in 



Ref. [24 



One of the goals for this study is the search of phenomena of the step dynamics in order 
that the mechanism of the crystallizing matter transfer should be determined. The distinct 
qualitative difference takes place for step bunching. At the surface diffusion in the presence 
of the heating current, step bunching appears both for step up and down direction of the 
current, whereas at the volume diffusion in the solution flow, the only steps moving parallel 
to the flow can be unstable. 

The lateral instability considered here occurs both at the surface diffusion due to the 
Schwoebel's effect and at the volume diffusion due to the drift of the diffusion field by the 
solution flow. However, in the case of the volume diffusion, the instability is controlled by 
the intensity of the solution flow and the flow direction with respect to the step motion. 
While interpreting experimental results, one should bear in mind the possibility that a 
lateral instability can be excited by the solution flow at the surface diffusion, because the 
delivery of the crystallizing units to the crystal surface goes through the moving solution at 
any rate. In the presence of the surface diffusion with the diffusion length \2d = (D2T S ) 1 ^ 2 , 
where D 2 is the surface diffusivity and t s is the characteristic time of the desorption of the 
units back into the solution, the whole path of the units includes diffusing in the solution to 
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a sink strip with the width of order of \2d and then diffusing on the crystal surface within 
this strip to the step. The drift of the volume diffusion layer by the solution flow leads to 
a nonuniformity of the supersaturation along the disturbed step, so the lateral instability 
is possible. In this case, only long wave perturbations, q < A/\2d, where A ~ 1, can be 
unstable, corresponding critical wavelength depending on neither the flow intensity nor the 
step velocity. Thus, these differences of the step dynamics may be the criterion determining 
the delivery mechanism. 

At the same time with the diffusion processes, the morphological instability can be con- 
ditioned by an impurity effect. The processes resulting in the loss of vicinal-surface stability 
of the prismatic faces of ADP and KDP under the action of the supersaturation nonunifor- 
mity and the presence of growth-decelerating impurities had been observed by in situ laser 
interferometry J29|. It had been established that the instability are most pronounced in the 
supersaturation range corresponding to the maximum partial derivatives of the step velocity 
with respect to the supersaturation and the impurity concentration. As it was found in Ref. 
an impurity influence can lead to the S-shape dependence of the step velocity on the 



supersaturation and on the step curvature too. At the supersaturation near by this hys- 
teretic range, the step is unstable against the lateral perturbations. A model of instability 
of the step spacings owing to impurity effect was considered in Ref. |yj . These instabilities 
conditioned by the impurities can be separated through the dependence on the impurity 
concentration. 
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FIGURE CAPTIONS 

Fig. 1. Scheme of singular crystal surface with single step. 

Fig. 2. Integation path in expression (]17D . 

Fig. 3. Scheme of the concentration isoline near by rectilinear step in the solution flow. 
Fig. 4. Dependence of Iq on dimensionless step height H = h/A. 
Fig. 5. Dependence of a{Q). 

Fig. 6. Dependence of / on dimensionless wave vector of lateral perturbation Q = qA. 

Fig. 7. Shaded zone corresponds to unstable lateral perturbation during contrary motion of 
step and solution flow. 
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